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We investigate the transport characteristics of a redox system weakly coupled to leads in the 
Coulomb blockade regime. The redox system comprises a donor and acceptor separated by an in- 
sulating bridge in a solution. It is modeled by a two-site extended Hubbard model which includes 
on-site and inter-site Coulomb interactions and the coupling to a bosonic bath. The current voltage 
characteristics is calculated at high temperatures using a rate equation approach. For high voltages 
exceeding the Coulomb repulsion at the donor site the calculated transport characteristics exhibit 
pronounced deviations from the behavior expected from single-electron transport. Depending on 
the relative sizes of the effective on-site and inter-site Coulomb interactions on one side and the 
reorganization energy on the other side we find negative differential resistance or current enhance- 
ment. 

PACS numbers: 71.27.a, 34.70.e, 82.39.Jn 



Electron transfer is fundamental in chemistry and bi- 
ology. In a polar environment the redox process in- 
volves strong coupling to the underlying nuclear motion 
and is usually dominated by the nuclear reorganization 
that accompanies the charge rearrangement. Such a re- 
dox system contacted to leads can serve as a component 
of molecular electronic devices^. Recent experiments^ 
showed that the electron transport in redox systems de- 
pends strongly on the environment (solvent) and can 
be described by a thermally activated electron transfer 
process mediated by thermal fluctuations of the environ- 
ment. 

Although the standard description focuses on single- 
electron transfer 3 multiple electron transfer is the dom- 
inant process in some chemical systems^^. Many- 
particle states have to be taken into account and the rate 
between these states depends strongly on the Coulomb 



interaction^. Contacting such a redox-system to leads 
and applying a voltage larger than the effective on-site 
Coulomb repulsion, e.g., on the donor, it can be occu- 
pied by multiple electrons. Then the single-electron pic- 
ture breaks down and the transport is affected by elec- 
tron correlations^. In the present paper, we study the 
current through such a redox sytem taking the solvent 
explicitly into account for the first time. The central 
results are displayed in Fig|2] and Figj3] where negative 
differential resistance (NDR) or current enhancement is 
be found, respectively. These effects occur at voltages 
when multiple electrons populate the DA system. 

We consider a two-site model where the redox sys- 
tem comprises donor (D) and acceptor (A) which are 
separated by an insulating bridge (see Fig. [T|). The 
system is embedded in a solvent and weakly coupled 
to metallic electrodes. The Hamiltonian is given by 
H = i^DA + Hda-c + Hq where the DA system is de- 
scribed by 
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FIG. 1: Schematic view of the redox system with donor (D) 
and acceptor (A) coupled to the leads L and R. The elec- 
tronic degrees of freedom of the DA system are coupled to 
the environment comprising internal vibrations and the sol- 
vent dynamics. The current is calculated as a function of the 
bias voltage V& and gate voltage V g . 



As written Ci a and q a denote annihilation and creation 
operators for fermions with spin a in an orbital local- 
ized at site i and = cl a Ci a is the number-operator. 
The first four terms correspond to an extended two-site 
Hubbard model, with on-site energies hopping matrix 
element A, on-site Coulomb repulsion U, and inter-site 
Coulomb repulsion Uda- The two last parts describe the 
electron-boson coupling (in the standard polaron form) 
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FIG. 2: (color online) Current I as a function of V& for 
C/eff = leV^ and [/da = 2eV (upper panel, black solid line) in 
comparison to the case of independent electrons (upper panel, 
red dashed line) as well as a contour plot of dl/dV vs. V g and 
Vb for correlated electrons (lower panel). The other parame- 
ters see in the text. Inset: Potential surfaces for the states: 
|TI,0), |T, I) and |0,TI>. 



and the free bosonic bath with boson creation and anni- 
hilation operators b^ n and b n , respectively. Furthermore, 
the coupling between leads and DA system is given by 
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while the dynamics of the free fermions in the leads is 
accounted for by H c = ^2 ky a;i=L y R e k n hk,a- In the fol- 
lowing, we focus on the Coulomb blockade regime with 
weak coupling between leads and DA system which is 
strongly coupled to the bosonic bath. The electrons tun- 
nel sequentially on and off as well as migrate via ther- 
mal activation through the DA system. At high tem- 
peratures, exceeding the activation energy, Marcus ap- 
proximation is valid where only one bosonic mode with 
frequency uoq is considered 7 . We obtain an effective 
model where the parameters are renormalized according 
to: cd, a = £d,a — E\, [/ e ff = U — AE\ (with the reorga- 
nization energy E\ = \ 2 /ujo, A — > A and V k — > Vfc. 

The DA-system has 16 basis-states \D,A) describing 
the occupation on D (A), with different number n = + 
via of electrons where ud(tia) is the number of electrons 
on D (A): | |,0>,| 1,0), |0,|>, |0,|>,| U>, I U>, I TI 
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FIG. 3: (color online) Current I as a function of Vb for 
Ueft = 3eV and V = leV (upper panel, black solid line) in 
comparison to the case of independent electrons (upper panel, 
red dashed line) as well as a contour plot of dl/dV vs. V g and 
Vb for correlated electrons (lower panel). The other parame- 
ters see in the text. Inset: Potential surfaces for the states: 
I t> is 0) 5 I T> 1) an d |0, ||). Note that in this case the transfer 
from | ||, 0) to | |,|) is in the activation- less regime. 



,o),io,u>, 1 iu>, 1 iu>, 1 r,u>,i uo,i u,u> and 

|0,0). 

Assuming fast dephasing and incoherent transfer the 
population dynamics can be described by kinetic equa- 
tions: 



P s — ^ k r ^ s P r ^ ^ k s ^ r P s , 



(5) 



where P s is the probability that the system is in state s 
and k r ^ s is the rate of the transition from state r to state 
s which is specified in the following. If the occupation ud 
is changed by one at fixed ua the rate is 

For the reverse process we find 
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(n D ,nA)^(n D — l,nA) 
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f (1 - f(AE sr - hl)) \{s\dr,\r)\ .(7) 



The hybridization strength f = 27r^ fe |T4| 2 ^(^ — e&) is 
assumed to be energy independent and equal for R and 
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L, / is the Fermi distribution and p is the density of 
states of the leads. We assume k B T ^> hT. Similar 
rates are obtained if ua is changed and is constant 
by exchanging L with R and D with A. The chemical 
potential of the lead L and R is defined with \±l,r — 
€f ± eVb, where Vb is the bias voltage. A finite gate 
voltage leads to Id — > £d + and — > + V^. 

The transfer rates, at conserved particle number n, is 
given by the Marcus theory in the nonadiabatic limit 
(A < u; ) 



k (nn,n A )^(n D ±l,n AT l) = 27T A 2 F( AE sr ) , 
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where the Franck- Condon integral 

1 



F(AE sr ), 



F(AE sr ) 
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depends on the energy difference AE sr between state 
s and r and the corresponding reorganization energy 
e.g., £?jT.»M0,T> = ^ ^ITi.oMi.t) = 4 ^ and 

E \n,o)~\o,in = mEx 

We calculate the stationary current (I) as a function 
of bias voltage Vb and gate voltage as 

i/io = Y.^7Z D+1 Pr - K^rr~ lp s) ( n ) 



exponentially on C/ e ff — Uda - ^E\- In Fig |2] and Figj3] 
(upper panel) we display / for /7 e ff = leV^UoA — 2eF 
and ZTeff = 3eV, Uda — leV, respectively, in comparison 
to the current where electron correlations are neglected 
(red dotted line) for V g = 0. For eV& < /7 e ff both curves 
are identical but at eVJ, > /7 e ff two electrons occupy D. 
The electrons are transferred via a transition from state 
I T1,0) to I T,|), | |,|) or |0,||) . The potential surfaces 
for these many-particle states are displayed in the inset. 
In Fig 12] the transfer is in the normal regime and in Fig. [3] 
in the activation- less regime since Z7 e ff — V = ^E\. In the 
latter case the current is enhanced for eV& > U e ^ = 3eV 
by a factor of 10000 in comparison to the uncorrelated 
case. For (U e s - V - 4,E X ) 2 > 4E x k B T the transfer is 
in the normal or inverted region and negative differen- 
tial resistance (NDR) is present as seen in Figj2] (lower 
panel). For our parameters (f7 e ff = leV, Uda = 2eV) the 
current is surpressed by a factor of 20. Furthermore, we 
examined a small dependency of the on-site energies on 
Vb and found that the results do not change qualitatively. 

In summary, we considered a solvated redox system 
weakly coupled to leads at high temperatures which is 
able to accept more than one electron. If the bias voltage 
is larger than the effective Coulomb repulsion, the donor 
can be occupied by two electrons and the transfer rates 
are dominated by correlation effects. We found that the 
current strongly depends on the difference between the 
on-site and inter-site Coulomb repulsion as well as the 
reorganization energy. The current is enhanced in the 
activation- less regime when the difference is zero, other- 
wise NDR behavior is obtained. 



by assuming P s = 0, P s = 1 and for the parameters 
l/(k B T) = 40eV, f = O.ObeV, A = O.leV, E x = 0.5eV. 
The current depends on the rate of electron transfer from 
D to A. Neglecting electron correlations two electrons 
on D are transferred each with a single-particle Marcus 
rate. However, electron correlations have to be taken 
into account (even for U = Uda = 0) for a correct de- 
scription, e.g., the rate for | ||,0) — > | |,|) depends 
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